Experimental Demonstration of Nonuniform Frequency Distributions of Granular Packings 
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We developed a novel experimental technique to generate mechanically stable (MS) packings of frictionless 
granular disks. We performed a series of coordinated experiments and numerical simulations to enumerate the 
MS packings in small 2D systems composed of bidisperse disks. We find that frictionless MS packings occur 
as discrete, well-separated points in configuration space and obtain excellent quantitative agreement between 
MS packings generated in experiments and simulations. In addition, we observe that MS packing probabilities 
can vary by many orders of magnitude and are robust with respect to the packing-generation procedure. These 
results suggest that the most frequent MS packings may dominate the structural and mechanical properties of 
granular systems. We argue that these results for small systems represent a crucial first-step in constructing a 
statistical description for large granular systems from the 'bottom-up'. 
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PACS numbers: 61.43.-j, 81.05.Kf, 63.50.Lm, 83.80.Fg 

The power of equilibrium statistical mechanics is that it en- 
ables the evaluation of macroscopic state variables (such as 
temperature and pressure) of a macroscopic system in ther- 
mal equilibrium simply by counting microstates. A number 
of recent studies have applied similar statistical methods to 
describe dense granular materials (H 0, 01. For example, 
Edwards-ensemble descriptions are based on an assumption 
that all mechanically stable (or 'jammed') configurations of a 
granular system under a given set of macroscopic constraints 
are equally likely. The Edwards' theory also implies the exis- 
tence of a temperature-like variable — the compactivity % fl. 

Despite the fact that granular media are dissipative and re- 
quire external driving forces (not thermal fluctuations) to ex- 
plore configuration space, there has been surprising success 
in describing these materials using statistical methods based 
on the Edwards' assumption ]5|]. For example, simulations 
of slowly sheared granular systems have shown that % can be 
used to quantify effective 'thermal' equilibrium, since differ- 
ent tracer particles achieve the same % Jo] and several equilib- 
rium measures of temperature all agree Q7|] . However, statisti- 
cal mechanics approaches for dense granular materials have 
been applied without directly testing the underlying funda- 
mental assumptions. In particular, the assumption of equal 
micrsostate probability has not been tested explicitly, and the 
relevant microstates have not been clearly defined. We ad- 
vocate a novel 'bottom-up' approach to constructing statisti- 
cal mechanics descriptions of dense granular materials — one 
where we enumerate the microstates and accurately measure 
the probabilities with which they occur. 

To do this, we performed a coordinated set of experimen- 
tal and computational studies of mechanically stable (MS) 
packings in small 2D granular systems undergoing vertical 
vibrations. To enable enumeration of MS packings (i.e., mi- 
crostates), we focused on small systems with no frictional 
forces. We show that in the absence of frictional forces, the 
set of MS packings is discrete; thus packing probabilities can 
be directly evaluated by counting the frequency with which 
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FIG. 1: (a) Schematic of experiment, (b) Coordinates (x c ,y c ) of the 
centroids of several N = 7 MS packings from experiments on plastic 
(squares) and simulations for Y p i ast i c (circles). The solid line shows 
the location of one of the centroids for 10~ 7 <J<Y*, where y* is the 
value at which the MS packing becomes unstable. The star indicates 
the location of the centroid at y* . (c) Probability distributions of the 
separation AR in configuration space between distinct MS packings 
(D) and between a given MS packing and the packing furthest away 
with the same contact network (C) for experiments (dashed lines) 
and simulations (solid lines), (d) Probability distribution for y* from 
simulations. The vertical line indicates y = Ypiastic- 

they occur in a long sequence of independent trials. In con- 
trast, frictional MS packings Hfl form continuous families, and 
therefore packing probabilities cannot be uniquely determined 
without first defining an appropriate probability measure. To 
generate frictionless packings in our experiments, we have de- 
veloped a novel technique where frictional forces are relaxed 
using small-amplitude, high-frequency vibrations. 

In both experiments and simulations we find the follow- 
ing four key results concerning the microstate distributions of 
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FIG. 2: (a) Number of distinct states N s found in N, trials for experi- 
ments on plastic (gray solid line) and N = 3 (dot dashed), 4 (dotted), 
5 (dashed), 6 (long dashed), and 7 (solid) simulations at y = Y p i as ti c • 
The horizontal (vertical) axis is scaled by the total number of MS 
packings A^ ot (trials, A', tot ) at saturation. The experimental curve was 
obtained by fitting _/V' ot and Af f tot to simulations, (b) Nf l vs. N from 
simulations (circles) and experiments (squares). The solid line has 
slope 1.2. (c) Sorted probability of MS packings for N = 7 vs. 
index (k + N l ° l — N s ) / N l ° l for simulations (solid) and experiments 
(dashed). N„/Nl ot ~ 0.695 for experiments was determined from (a). 

granular packings in the zero-friction limit: 1) There exist a 
finite number of discrete MS packings that grows exponen- 
tially with system size. 2) The frequency with which these 
packings occur is highly nonuniform. 3) The sets of packings 
found in experiments and simulations of model granular me- 
dia are very similar. 4) The packing frequency is relatively 
insensitive to the packing preparation protocol. We argue that 
the above important new results need to be incorporated into 
statistical descriptions of dense granular media. 

Experiments A schematic of the apparatus used to gen- 
erate MS frictionless disk packings is shown in Fig. \V[a). 
A mixture of thin disks of thickness 3. 175 ±0.003 mm and 
two different diameters a s and O/ were confined between two 
glass plates separated by 3.20 ± 0.01 mm and rested on a thin 
plunger connected to an electromagnetic shaker through a slot 
in the bottom of the cell. The shaker enabled us to apply ver- 
tical vibrations at variable amplitude and frequency to repeat- 
edly generate static particle packings. 



The particle mixtures consisted of (N+l)/2 small and 
(N — l)/2 large particles, with N = 5 and 7. For these sys- 
tems we enumerated the majority of frictionless MS packings 
by performing A?, = 1.8 x 10 5 (N = 7) and 1.2 x 10 4 (N = 5) 
independent trials using the protocol described below. The ra- 
tio of disk diameters was d = 0//a s = 1.2520 ±0.0003, and 
the ratio of the cell width L to the small particle diameter was 
X = L/a s = 4.253 14 ±0.00001 (2.65 ±0.02) for N = 7 (5). 
We used bidisperse systems to prevent ordering. 

The degree to which our particles behave as hard disks 
can be estimated using the dimensionless stiffness parameter 
y = m s g/ka s , where g is the gravitational acceleration, m s is 
the mass of a small particle, and k is the effective spring con- 
stant of the elastic interparticle interaction. By measuring the 
deformation of single plastic (steel) disks under gravity, we 
estimate y p i as tic ~ 1.85 x 10~ 3 (y ste ei ~ 3 x 10~ 7 ), which im- 
plies that the deviation of the particle packings from hard-disk 
behavior is small (cf., Figs. |TJi>) andQJcO). 

To generate an ensemble of frictionless MS packings, we 
repeatedly performed the following protocol: The plunger 
was first oscillated at high amplitude and low frequency (50 
Hz) for 100 ms to randomize particle positions. The system 
was allowed to relax under gravity with the shaker turned off 
for 400 ms. We then applied a low-amplitude, high-frequency 
(400 Hz) oscillation for 500 ms, which excites particle rota- 
tion and relaxes frictional particle-particle and particle-wall 
interactions. Finally, the oscillations were turned off and po- 
sitions of particle centers were determined to an accuracy of 
A/o v = 6 x 10~ 6 using a digital camera and particle-tracking 
software. The MS packings in experiments were identified us- 
ing the set of particle positions Rj = {r\ ,7i, . . . for each 
configuration i, where r, are the x- and y-coordinates of the N 
particles. (See Fig.[T]a).) 

Computer Simulations We also performed molecular dy- 
namics simulations of gravitational deposition of bidisperse 
frictionless disks. Our goal was to determine how key features 
of MS packing distributions depend on the particle-deposition 
process, and to identify which features are robust, i.e., do not 
depend on specific details of the dynamics. Thus, in the simu- 
lations we do not exactly mimic the packing-preparation pro- 
cess in experiments. In particular, we do not model frictional 
contact forces, but instead we use velocity-dependent resis- 
tance forces to dissipate energy. However, geometrically sim- 
ilar sets of MS packings are needed for a detailed comparison 
of the probability distributions. Thus, we closely match the 
cell width X, particle size distribution, gravitational force, and 
elastic interactions in simulations and experiments to obtain 
very similar sets of MS packings. 

We assume that the disks interact via a finite-range, purely 
repulsive linear spring force 

F r (r i j) = —8 i j&(-8 i j)? ij , (1) 

which mimics elastic interparticle repulsion. £ is the char- 
acteristic energy scale, ry is the separation between particles i 
and j, Ojj = (Ci + Cj) /2 is the average diameter, 8y = r,y — an 
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is the interparticle overlap, ry is the unit vector connecting 
particle centers, and &(x) is the Heaviside step function. 

To create each MS packing, we randomly place particles in 
a square cell of size L, with no particle overlaps. We then al- 
low initially stationary particles, interacting via elastic repul- 
sive forces (Q~|i and dissipative forces proportional to relative 
particle velocities, to fall under gravity. The system evolves 
according to Newton's equations of motion 



N 

mi3i = -nijgy + Y, F r (nj) - bQ(-8ij)vij ■ hj 

Mi 



fij+Fr, (2) 



where a\ is the acceleration of particle i, v,- ; - is the relative 
velocity of particles ; and j, and b is the damping coeffi- 
cient. The particle-wall interaction force F w has an analogous 
form to the particle-particle interaction (0, with energy scale 
E w = 2s. We set the dimensionless damping coefficient to 
b = a s b/y/m s E = 0.25. The simulations are terminated when 
the total force F tot on each particle is vanishingly small. (In 



most simulations we used the threshold F tat < F m 
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To determine which of the relaxed configurations are me- 
chanically stable, we calculated the eigenvalues of the dynam- 
ical matrix @|. The MS packings possess 2N' positive eigen- 
values, where N' = N — N r , and N r is the number of 'rattler' 
particles. Rattlers have fewer than three contacts (including 
wall contacts) and if present typically rest on the bottom. In 
the simulations, we distinguish distinct MS packings by com- 
paring the eigenvalue lists. The eigenvalues are considered to 
be equal if they differ by less than the noise threshold 10~ 6 . 
Less than 1 % of the distinct MS packings contain rattlers. In 
these configurations, we ignore the translational degeneracy 
of the rattlers — two configurations with the same contact net- 
works of non-rattler particles are treated as the same. To com- 
pare simulation and experimental data in configuration space, 
we omit the rattler particles, and consider only positions of 
particles forming the contact network. To enumerate all MS 
packings and accurately measure their frequencies, we consid- 
ered small systems in the range N — 2tol particles. Systems 
with an even number of particles contained equal numbers of 
large and small particles, while systems with an odd number 
contained one additional small particle. 

Results We find that MS packings in frictionless granular 
systems occur as discrete points in configuration space Il4ll . 
In Fig. \V[b) we display the coordinates (x c ,y c ) of the cen- 
troids of MS packings in a small region containing several 
microstates for N = 7 (plastic discs in experiments and par- 
ticles with Ypiastic in simulations). The results show that the 
MS packing centroids are indeed distinct and well-separated. 
Moreover, the experimental and simulation points agree. Our 
simulations do not involve static frictional forces, thus this 
agreement indicates that our novel experimental technique is 
able to generate mechanically stable frictionless packings. 

Fig 02c) shows that the experimental and numerical scat- 
ter in the MS packing centroids is several orders of magni- 
tude smaller than the average separation between discrete MS 
packings in configuration space. We find that the average dis- 




FIG. 3: (a) Distance AR„ and AR„„ between each MS packing from 
experiments on plastic disks and the nearest (dotted) and next-nearest 
(solid) MS packing from simulations at Ypiastic vs - index k for the ex- 
perimental MS packings sorted in order of increasing AR„. The ver- 
tical line at k = 618 separates matched MS packings to the left from 
those that are unmatched, (b) The probability with which MS pack- 
ings occurred in experiments, /^ xp versus probability of the matching 
state in simulations P?™. 



The solid line has slope one. 



tance between distinct MS packings in configuration space 



is approximately 10 G s , whereas the maximum size of the 
scatter is 10~ 12 G S in simulations and lO^CJj in experiments. 
Based on this observation, in our analysis of experimental 
data, two packings Rj and Rj are considered to be the same 
microstate if AR = \Rj-Rj\/o s < AR* = 0.01. 

To determine if particles in our system can be treated as 
hard disks, we tested the sensitivity of the numerically gen- 
erated MS packings to variation in the stiffness parameter y. 
The solid line in Fig.QJ^) shows the change of position of the 
centroid of a MS packing when y is increased from 10~ 7 to the 
critical value y* where a sudden change in the particle contacts 
occurs (y* =0.06 for this particular packing). While the over- 
all variation of the position of the centroid is significant, the 
position of the centroid for y = Ypiastic is essentially indistin- 
guishable from the position in the hard-disk limit y — > 0. The 
distribution of the critical parameters Y* for the set of all simu- 
lation MS packings is depicted in[TJ^)- These results indicate 
that most of the hard sphere packings remain stable even when 
the stiffness parameter is increased above Ypiastic- 

In simulations, we are able to perform an extremely large 
number of trials and find nearly all MS packings in small sys- 
tems. In Fig. rfla), we show the number of distinct states 
N s as a function of the number of trials N t for systems with 
N = 3 through 7 particles. In all cases (except = 7), we 
saturate the packing-generation process in the sense that we 
do not generate new MS packings when the number of trials 
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is increased by a factor of 10 beyond A^ tot . In experiments 
we followed a similar procedure, but we were unable to fully 
saturate the curves because of insufficient number of trials. 
Fig. shows that in experiments and simulations the total 
number of MS packings grows exponentially with the system 
size, Nl° l ~ e . The exp onent a « 1.2 is the same as found 
for periodic systems 114J1 ; however, the prefactor is larger by 
roughly an order of magnitude. The number of trials required 
to reach saturation of the simulation packing-generation algo- 
rithm (cf., Fig. 0a)) also grows exponentially with N, but with 
a larger exponent. Since A',' ' grows rapidly with system size, 
enumeration of A^ tot ~ 728 MS packings for N = 1 requires 
~ 10 9 trials. This large number of trials stems from the 
extremely nonuniform packing probability distribution. 

The frequency distributions for MS packings are extremely 
nonuniform in both simulations and experiments. As depicted 
in Fig. |2 c) the MS packing probabilities vary by many orders 
of magnitude. In addition, we find quantitative agreement in 
the shape of the frequency distributions, which implies that 
frequencies of the MS packings are only weakly sensitive to 
the dynamics used to generate them. 

To make a quantitative comparison between MS packings 
found in experiments and simulations, we calculated the dis- 
tance in configuration space AR between each MS packing 
generated in experiments and the nearest and next nearest 
MS packings found in simulations. In Fig. 0a), we show 
the nearest-neighbor and next-nearest neighbor separations 
AR n and AR n „ for experiments with N = 1 plastic disks and 
the corresponding simulations at Yplastic versus index k sorted 
by increasing AR n . For approximately 80% of the packings 
the nearest-neighbor and next-nearest-neighbor distances are 
well separated, with AR„ much smaller than the average dis- 
tance between packings shown in Fig. []I C )- This separation 
of length scales in configuration space allows an unambigu- 
ous match between packings found in experiments and simu- 
lations. At k w k* =618 the distance AR„ rapidly increases, 
and packings with k > k* are unmatched. 

We find that n = 110 (m = 146) MS packings from simu- 
lations at Ypiastic (experiments on plastic disks) are unmatched 
for N = 1 . However, approximately 70% of the unmatched 
experimental packings are unstable when used as initial con- 
ditions in the simulations. (In experiments, these packings are 
likely stabilized by residual frictional forces or small inaccu- 
racies in the numerical representation of the experiment.) For 
example, the results shown in Fig. \\\d) indicate that stabil- 
ity of some packings may depend on the stiffness parameter 
near Ypi ast i c . The unmatched MS packings from simulation 
most likely result from insufficient experimental statistics. 

To determine the sensitivity of packing probability distribu- 
tions on the particle deposition process, in Fig. 0£>) we com- 
pare the experimental probabilities P^ xp with the probabilities 
found in simulations P| lm for the set of matched states. We 
demonstrate a strong correlation between P^ xp and P^ lm : likely 
packings in experiments tend to be likely in simulations, and 
rare packings in experiments tend to be rare in simulations 



(although there is also a significant scatter). In fact, we cal- 
culate that the rms deviation in the probabilities of matched 
packings in simulations and experiments is only 8% of the 
probability of the most frequent MS packing. Since the dy- 
namics in experiments and simulations is quite different, this 
is an important result, which implies that properties of static 
frictionless packings are weakly dependent on the packing- 
preparation protocol. 

Conclusions We introduced a novel experimental method 
to generate frictionless MS packings of granular materials. 
This method is crucial for studies aimed at differentiating the 
effects of geometrical constraints and friction on the struc- 
tural and mechanical properties of jammed granular systems. 
We performed coordinated experimental and computational 
studies of frictionless MS packings in small systems, which 
showed that MS packing probabilities are extremely nonuni- 
form and relatively insensitive to the procedure used to pre- 
pare them. In future studies, we will investigate the conse- 
quences of our present results for the microstate statistics in 
macroscopic granular systems treated as a collection of nearly 
independent small subsystems. We will also dial in frictional 
contacts to generate continuous geometrical families of MS 
packings that occur even at fixed y, and then compare the 
statistics of these continuous sets of packings to that for dis- 
crete MS packings. 
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